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The  finite  temperature  generalization  of  the  Quinn, 
Adler  theory  for  the  decay  rate  of  electronic  excitations  in 
a  normal  Fermion  system  is  presented.  The  theory  is 
developed  in  terms  of  quasi-particle  parameters  and  without 
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the  dielectric  response  of  the  system.  The  relations  derived 
are  pertinent  to  the  lifetime  problem  in  real  solids;  in  par¬ 
ticular,  local  field  effects  are  rigorously  taken  into  account. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


CONTENTS 

Introduction  1 

General  Theory  2 

Electron  Self-Energy  5 

References  14 

APPENDIX  A  —  An  Approximate  Evaluation  of  Finite 
Temperature  Corrections  for  a 
Simple  Case  1  5 

APPENDIX  B  —  Spectral  Density  Functions  and 

Sum  Rules  19 


IV 


Finite  Temperature  Theory  for  the  Attenuation  of  Quasi- Par  tide 

Excitations  in  Real  Metals 

INTRODUCTION 

Several  years  ago  Adler  ^  presented  a  theory  for  the  attenuation 

rate  of  hot  electrons  in  metals.  The  theory  was  developed  along  lines 
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similar  to  the  earlier  work  of  Quinn  and  Ferrell,  and  Quinn  concerning 
the  properties  of  an  interacting  electron  gas.  All  of  the  above  calculations 
were  carried  out  for  the  case  of  zero  temperature.  As  far  as  we  know, 
the  corresponding  extension  of  the  theory  to  finite  temperature  has  not 
been  given  in  a  general  context.  It  is  the  purpose  of  this  note  to  simply 
provide  such  a  generalization,  in  a  form  relevant  for  band  electrons  in 
solids.  Let  us  comment  at  once  that,  for  the  case  of  truly  hot  electrons, 
finite  temperature  corrections  cannot  be  very  important.  It  is  only  for 
the  case  of  low-lying  excitations  that  these  effects  need  be  taken  into 
account.  Furthermore,  because  we  are  concerned  only  with  low  energy 
excitations,  we  may  consistently  formulate  the  theory  in  terms  of  quasi¬ 
particle  energies. 

For  the  most  part,  the  results  presented  in  this  note  are  formal. 

For  the  case  of  free  electrons  at  sufficiently  low  temperature,  it  is 
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known  that  the  decay  rate  for  low-lying  quasi-particle  excitations  must 
have  exactly  the  form 

rk(£)  =  CkC(£^)2  +  Tr2(kBT)2]  * 

A  brief  discussion  of  this  result,  within  the  framework  of  a  simple  RPA 
treatment  is  presented  in  Appendix  A.  Qualitatively,  one  expects  to  find 
a  similar  result  for  the  case  of  normal  simple  metals.  In  particular,  the 
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results  for  low  energy  excitations  obtained  by  Adler,  barring  the  more 
pathological  cases,  e.  g.  ,  a  cylindrical  Fermi  surface,  which  he  has 
discussed,  may,  with  suitable  assumptions,  be  shown  to  generalize  at 
finite  temperature  in  a  manner  analogous  to  that  presented  in  Appendix  A. 

Finally,  in  real  solids,  one  expects  the  damping  due  to  electron- 
phonon  interaction  to  be  comparable  in  magnitude  to  the  effects  we  are 
considering  here.  This  problem  may,  however,  be  handled  in  a  parallel 
manner.  A  rather  thorough  discussion  of  the  electron-phonon  interaction 
at  finite  temperature,  based  on  a  model  of  free  electrons  interacting  with 
longitudinal  acoustic  phonons,  may  be  found  in  the  book  by  Abrikosov,  et 


GENERAL  THEORY 

Our  theoretical  approach  is  based  on  the  well-known  finite  temper- 

ature  perturbation  theory  as  developed  by  Luttinger  and  Ward,  ^  and 
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Luttinger.  ’  In  this  section  we  give  a  brief  review  of  the  relevant 
propagator  formalism. 

We  take  the  Hamiltonian  to  be  of  the  standard  form 


H  =  H  +  H '  , 
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with 
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where  v(q)  =  4ne  /q  and  p(q)  >  -  (r|e  |r/)  . 

r  r  1 

Next  we  define  the  exact  single  particle  propagator  and  polarization 

f  9,  10  f.  ’ 

propagator  respectively  as 


g'_/(p',p')  =  Tr{e'e{H‘MN"n)T[a  W) „a  (p')H]}  ,  (4) 

r  xi  r  ri 
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where  P  ^  3  >  P  ^  0  and  all  operators  are  in  the  Heisenberg  representation, 
e.g- , 


p*H  -g'H 
p5/(e)H  =  e  p3,e 


■r 


rr 
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(6) 


Making  use  of  the  cyclic  invariance  of  the  trace  in  these  expressions 
one  obtains  the  Fourier  representations 


g  /(8'>  p")  =  7  y  s'  /(z  )e 
rr  p  L  rr  -C 


«4(p'-p#) 


(7) 


p—  /G3',e")  = 
qq 
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where  z,  =  tt i (2 -C  +  1 ) / 3  +  p,  §.  =  Z Tri-C / 3  and  takes  on  all  integer  values. 

In  perturbation  theory,  S/  ,  ( z  )  is  given  by  the  sum  of  all  electron 

r  r  '0 

propagator  diagrams  in  which  a  Fermion  line  labeled  r'  enters  and  a  line 
labeled  r  emerges.  Likewise,  II )-*-♦,  is  given  by  the  sum  of  all 

i  qq 
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—  t  f 

polarisation  diagrams  in  which  an  interaction  line  labeled  q  enters  and 
a  line  labeled  q  emerges  (the  external  lines  not  being  part  of  II7).  Further¬ 
more,  S  and  II/  satisfy  their  respective  Dyson  equations 


S  ,(z  )  =  s  (z  )  6  ,  +  s  (z  ) 

rr  u  r  -C  r  r  r  L* 
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where  G  is  the  proper  self-energy  function  and  II  is  the  proper  polarization 
function.  Now  we  may  label  the  states  [  r)  as  |kYc)  where  a  is  a  spin 
index,  Visa  band  index,  and  where  the  wavevector  k  is  confined  to  the 
first  zone.  In  this  reduced  zone  scheme,  S7  and  G  will  be  diagonal  in 
the  wavevector  index,  but  will  have  off-diagonal  elements  in  the  band 
indices.  Likewise,  from  quite  general  symmetry  considerations,  II 7  and 
11  have  off-diagonal  elements  only  for  q*s  differing  by  a  reciprocal  lattice 
vector. 


Finally,  we  may  write  down  the  very  useful  spectral  represen¬ 
tations  of  the  propagator  s  S  7  ,  (z  .  )  and  II  7  (Z.  which  are 

rr  -l  t  q  q 
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rr  -C  J 


dE 


A<E>rr' 

VE 


(11) 


cc. 

n'V55'=J 


B(E )-«-«/ 

dE 


(12) 


These  relations  are  obtained  in  the  usual  manner  by  evaluating  (4)  and 
(5)  in  the  representation  |n)  of  exact  eigenstates  of  H,  and  then  inverting 
(7)  and  (8)  to  obtain  the  corresponding  Fourier  coefficients.  The  spectral 
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density  functions  A  (E)^  ,  and  B(E)->^,  are  given  in  Appendix  B,  along 
with  their  various  general  sum  rules.  One  important  property  of  the 
spectral  density  functions  is 

A(E)  *  =  A(E)  #  ,  (13) 

rr  r  r 

j»{ 

B(E)-*-4/  =  B(E)->,-.  ,  (14) 

q  q  q  q 

i.e.  ,  the  diagonal  elements  are  real. 

ELECTRON  SELF -ENERGY 

We.  now  proceed  to  calculate  the  electron  self-energy  in  the 
approximation  which  retains  only  the  lowest  order  vertex  function.  In 
this  approximation,  one  obtains,  using  the  well-known  graphical  rules, 

G  ,(z.)  =  G(1),  +  G(2),(z.)  ,  (15) 

rr  v  rr  rr  v 


,  rd) 

where  u  / 
rr 


is  a  constant,  independent  of  z  ,  and  is  given  by 


[(rs|vjv's')  -  (rs|v  |s'r')] 
s  s 


1  Z  Ss's  (V,e 


V°+ 


t' 


(16) 


This  term  is  clearly  a  generalization  of  the  usual  Hartree-Fock  self¬ 
energy,  and  is  intimately  connected  with  one  of  the  sum  rules  satisfied 
by  the. electron  spectral  density  function  (see  Appendix  B).  The  second 
term  in  Eq.  (15)  is  more  interesting.  This  term  is  usually  associated 
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with  electron  correlation,  and  is  given  by 


v(ci)v(4')^  p_ /(q) p„_/(-q') 


r  s 


sr 


qq 


ss 
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Now,  as  Luttinger  has  discussed,  we  should  seek  a  similarity 

transformation  which  diagonalizes  the  matrix  C  6  /  +  G  ,(z0).  This 

r  rr  rr  t 

is  almost  surely  unfeasible  to  obtain  in  practice,  and  perhaps  even  in 
principle,  i.e.  ,  it  is  not  known  whether  such  a  transformation  even 
exists  in  general.  Therefore,  in  order  to  obtain  something  tractable, 
we  shall  proceed  by  simply  assuming  that  the  off-diagonal  elements  of 
G  (and  hence  also  of  S  / )  may  be  neglected. 

Because  G  ^  is  real,  this  term  does  not  contribute  directly  to 

r 

the  electronic  attenuation  rate,  i.e.  ,  its  effect  is  simply  to  modify  the 

quasi-particle  dispersion  spectrum.  Turning  to  the  consideration  of 

(2) 

G  (z  ),  we  may  invoke  the  spectral  representations  of  the  propagators 
r  v 

S'  and  n7  and  obtain 
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The  5.  /  sum  in  this  expression  is  now  trivially  done  using  the 
standard  prescription 
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(19) 
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where  N(§)  =  (e  -1)  is  the  Bose  function,  and  F(^)  has  only  simple 
poles.  Using  this  prescription,  we  have 


i  y  _j _ i _ 

*  4  V'E  V_VE 

_  f4(E)  N  '  (E  ' )  +  f  ~  (E)  N(E  ' ) 
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with  N'(E')  =  N(E')  +  1.  Having  carried  out  the  §  /  sum  we  may  now 

(2)  ^ 

continue  G  (z .  )  off  the  set  of  points  z,.  The  (unique)  continuation 

which  has  no  essential  singularity  at  infinity  may  then  be  written 


C^Z\z)  =  ~  V  v(q)v(q')7  p  (q)p  (-q') 
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We  note  that  the  function  defined  by  (21)  is  analytic  except  on  the 

(2)  *  (2)  * 

real  z-axis.  Furthermore,  from  (13)  and  (14),  we  see  G  (z)  =  G  (z  ), 

r  r 

i.e.  ,  there  is  a  pure  imaginary  discontinuity  across  the  real  z-axis.  Thus, 
for  z  near  the  real  axis,  we  have 
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and  hcncc 
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From  the  spectral  representation  of  II  7(§  ,)-<-*  ,  ,  it  is  then  easy  to  see 

l  qq 

that  this  result  may  be  written 
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The  functions  G  (eTis)  correspond  to  the  well-known  advanced  and 


(24) 


retarded  self-energy. 
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Now  for  e  near  |i,  the  statistical  factors  in  (24)  constrain  E  to 
values  near  p,  and  for  such  values,  we  assume  that  A(E)_  =  I  ImS^(E-is) 
is  a  sharply  peaked  function,  i.e.,  we  assume 


ir 


A(E)  =■*  6[E-e  -A  (E)] 
s  s  s 
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where  A  is  the  real  part  of  the  self-energy, 
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and  E  is  the  (we  assume  unique)  quasi-particle  energy  given  by  the 
s 

solution  of 


E  -  e  -  A  (E  )  =  0  .  (26) 

s  s  s  s 

As  always,  the  ultimate  justification  of  Equation  (25)  must  be  the 
self-consistency  of  the  final  result. 

Using  this  approximation,  we  obtain  for  the  imaginary  part  of  G 
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Now  Im  G  (e  *fis)  and  Im  G  (e-is)  are  related  to  the  damping  rate  for 
r  r 

particle  and  hole  propagation  respectively.  In  either  case,  because 
these  functions  differ  only  in  sign,  we  shall  simply  assess  the  lifetime 
for  the  corresponding  excitation  through  the  relation 
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If  one  likes,  the  above  result  may  now  be  expressed  in  terms 
of  the  so-called  dielectric  response  function.  Adler  ^  has  discussed 
this  function  in  some  detail  within  the  R  PA  formalism.  Let  us  define 


K(e-is)-.-.,  =  6-*  +  II  (e -is  )-*-*,  v(q* )  ,  (30) 

qq  q»q  qq 


and  the  dielectric  response  function  K  ^(e-is)-»-»,  by  the  relation 

q  q 
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Then  from  analytic  continuation  of  (10),  we  obtain 
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(32) 
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Making  use  of  this  result  in  Eq.  (29),  we  find 


r  (E  )=Im2.  V  )'  Z  P  (q)p  (-q')v(q) 
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This  result  is  most  simply  written  out  explicitly  in  the  extended 
.one  scheme;  one  obtains 
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In  this  expression  we  have  suppressed  spin  indices,  K  and  K  are 
reciprocal  lattice  vectors.  The  matrix  element  M  in  (34)  is  given  by 

+>  — ♦ 

— »  — ►  — ►  p  3-4  \  K  •  r  — *  —4 

M£(k,k-q)=J  d  re  u^(r)  u£_-(r)  ,  (35) 

unit  ^ 

cell 


11 


where  the  u's  are  the  usual  periodic  part  of  the  Bloch  functions, 
normalized  to  unity  in  a  primitive  cell. 

At  zero  temperature,  (34)  reduces  to  give 


^k^k)  =  "V"  L>  Z  Zk-q  +  ^  k-q) 

q  KK' 
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for  the  case  >  E  .  and 
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for  the  case  E^  <  E  .  Remark  that,  aside  from  the  renormalization 
k  f 

factor  Z  and  the  appearance  of  quasi-particle  energies  in  place  of  the 
unperturbed  single -particle  energies,  relation  (36a)  is  exactly  Adler  *  s 
zero  temperature  result  for  the  decay  rate  of  a  particle  introduced  above 
and  Fermi  surface.  The  interpretation  of  (36b)  is  equally  clear,  namely 
it  represents  the  decay  rate  for  a  hole  created  below  the  Fermi  surface. 

To  conclude,  we  should  like  to  make  several  comments  concerning 
the  dielectric  response  function.  As  Adler  has  discussed,  the  off- 
diagonal  elements  of  the  dielectric  response  function  give  rise  to 
umklapp-local-field  corrections.  Originally,  on  the  basis  of  a  one-OPW 
band  approximation  for  aluminum,  and  using  the  RPA  dielectric  response, 
Adler  estimated  that  umklapp  and  umklapp-local-field  effects  might  lead 
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to  a  correction  of  30%  or  smaller  for  the  lifetime  of  low-lying  excitations, 

while  these  effects  might  be  quite  important  in  the  plasma  production 
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region.  Recently,  however,  he  has  retracted  these  estimates  and  now 
allows  that  these  effects  probably  make  a  rather  small  correction  in 
both  regions.  Thus,  except  perhaps  for  certain  special  cases,  it  would 
seem  that  there  is  justification  for  ignoring  corrections  arising  from  the 
off-diagonal  elements  of  the  dielectric  response  function.  In  this  case, 
we  may  simply  write 


K  *(e  -is)-»-, 

qq 


(37) 


K(q, G  -is) 


where 


K(q,e-is)  =  1  +  v(q)  II(c-is)-~» 

q  q 


(38) 


is  the  usual  longitudinal  dielectric  constant.  Using  tins  approximation, 
(34)  reduces  to  give 
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APPENDIX  A 


An  Approximate  Evaluation  of  Finite  Temperature  Corrections 

for  a  Simple  Case 

We  wish  to  give  an  approximate  evaluation  of  (34)  for  the  case  of 
free  interacting  electrons.  We  further  restrict  ourselves  to  the  RPA 
dielectric  response  of  the  electron  gas,  and  shall  simply  assume  that 
the  quasi-particle  spectrum  may  be  approximated  by  the  non-interacting 
dispersion  law.  At  zero  temperature,  the  result  we  obtain  reduces  to 
that  quoted  by  Quinn. 

For  the  case  of  free  electrons,  the  result  (34)  becomes  simply 


rE(eE>  =  “~3  J  d3S  [£+<ek-?>  +  N<ek-  ek-5):l 

(Ztt  ) 


v(q)  Lm 


1 


(A.  1) 


K(q, 


_  dS  — »  dG  7~* 

t  .  ,3-*  q  k-q 

Letting  d  q  =  — ^ J - ■ 


t-= - j —  ,  where  dS~*  is  an  element  of  surface  area  on 

IVe-31  q 


the  surface  €-»  -»  =  constant,  we  obtain 
k-q 


CD 

rk(ek»  ---hi  dek-Q  [£+(Ek-J>  +  N(Sk-  ek.?> 


(2ir)  o 


I 


dS-*v(q)  Im  K(q ,  -  €-q  is ) 


(A.  2) 


=  cons 

k-q 


t  T^Gk-ql  |K(q,  e-- £-_--is)| 
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Now  for  the  case  M*  °f  low  energy  excitations,  the  statistical 

factors  in  (A.  2)  insure  that  the  energy  argument  of  the  dielectric  constant 

is  small  in  the  region  where  the  integrand  is  non-negligible .  Then 

* 

using  some  properties  of  the  well-known  RPA  dielectric  constant,  we 


have 


k  k2 

f  s 


Im  K (q ,  e-is)  =  -  - ~  , 

i  q 


(A.  3) 


|K(q,  G -is)  |  2  -  K(q,0)2 


(A.  4) 


for  sufficiently  small  G  and  q.  In  the  expressions  above,  k  is  the 

s 

inverse  Thomas-Fermi  screening  length.  Making  use  of  these  results 
in  (A.  2)  we  obtain 


rk(ek> = 


2,  l2 

e  k.k 
f  s 

8ttG  „ 


i 


de-*  -  (e-*-  e-*  -)  [f+(e-*  -)  +  N(e-*-  e-  -)] 
k-q  k  k-q  k-q  k  k  -  q 


(A.  5) 


It  may  readily  be  seen  that  we  are  justified  in  using  the  zero  temperature 
dielectric  constant  in  this  calculation. 
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Upon  carrying  out  the  surface  integral,  we  find 


W  UrXT-Xfe-r 


UJ 

f  de-.  (£r--£r;  — )  [ f  (er;  ->  +  N(e— -  e-*  -)] 

J  k-q  k  k-q  k-q)  k  k-q 


(A.  6) 


[ 


-1  ,  ksq  -k  +  P 

tan  k  +7T~T 


1 


S  kg  +  q  |k-P  | 


where  P  =  (  ~~  e-*  .  The  remaining  integration  over  e_*  is 

V  f  Z  k-qy  k-q 

clearly  quite  complicated.  On  the  other  hand,  for  the  principal 

contribution  to  the  integral  comes  from  the  region  ~  |a,  and  for 

k  — q 

such  values,  the  last  factor  in  (A.  6)  is  slowly  varying;  hence  we 
approximate  (A.  6)  as 


e2k. 


\  (  f  ^ 

\i-£\ 

JL[ 

'  'vk  y 

'  V  k  ) 

2  L 

2k  2k  k 
f  +  f  S 


ks  k^  +  4  k  f 
s  f 


]l(G^)  ,  (A.  7) 


with 


UJ 

I(eE>  =  J  d£k-5  (Ek-  Ek-?>  [f+(£k-q>  +  N(Ck-  £k-5)] 


(A.  8) 


Then  I(e_*)  evaluates  to  give 

K 

I(eK)  =  T  +  ir2(kBT)2  +•••]»  (A.  9) 

where  the  terms  which  have  been  neglected  vanish  exponentially  as  T  -»  0. 
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Our  ^inal  result  for  €  ~  \l  then  becomes 

K, 


rk(ck)=  h  Ce2kf)C"ir') 


2k  k 
f  s 


k  _  .  2k 

Y-£h!  +  , 

s  ks  k2  + 4k2 


] 


•Ctt)-T  C(ce-m)2  +  ,2(kBT)2]  . 

e 


2  .  2 


(A.  .10} 


In  a  similar  manner,  by  inserting  a  factor  of  —  (e-*  -  e-*  —  )  into 

n  k  K— c[ 

the  integrand  of  (A.  1),  we  may  obtain  an  approximate  expression  for 

13 

the  rate  of  energy  loss  for  an  incident  particle.  The  calculation 
proceeds  in  complete  analogy  to  the  one  above;  we  find 


de- 


dt 


•  is  ‘XX’CirX1*1’ 


-lfkf  +  2kf  ks 
ks  k2  +  4k.2 


] 


[(€r^)3  +1,2(kBT)2(ej-M)]  , 

e 


(A.  11) 


so  that  again,  the  lowest  order  finite  temperature  corrections  are  of 


order  T  . 
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APPENDIX  B 


Spectral  Density  Functions  and  Sum  Rules 

Making  use  of  relations  (4)  through  (8),  and  the  definitions  (11)  and 
(12),  it  is  straightforward  to  obtain  the  spectral  density  functions  for 
the  exact  one -electron,  and  polarization  propagators;  they  are,  respectively 


A(E)  , 

rr 


=  [eP(E‘M)  +  1]T  e"P(En"^Nn_n)  (nlah  |m)<m|a  |n)6[E-(E  -E  )]  , 

1  r  1  1  r 1  n  m 

(B.  1) 


nm 


B(E)-.-,, 

qq 


=  [e^E-  1]  L  )  e  ^En"uNn"^  (n|  p-,  |  m)  <m|  p  -|n>  6[E-(E  -E  )]  . 

v  Lj  q  -q  n  m 

(B.  2) 


nm 


It  may  readily  be  verified  that  B(E)--*,  satisfies  the  following  two  sum 


rule: 


u 

J 


dEB(E)-*-* ,  =-t  Tr  p  p-,])  =  0  , 

qq  v  -q  q 


(B.  3) 


J  dEEB(E)35,  =  -■i-Tr{e-e(H-wN-CI)[p_-[p3„H]]} 


q,  q 


(B.  4) 


—  (q.  q  ) Y  Trie  p^-}  . 


— ♦  “+  f 

W e  nots  that  the  right-hand  side  of  (B.  4)  vanishes  unless  q  and  q  differ 
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by  a  reciprocal  lattice  vector,  and  for  the  case  K  =  0,  we  obtain  the 
well-known  f- sum  rule 


I 


dEEB(E)- 


.2  2" 

n  q  n 


qq  m 


(B.5) 


Together  (B.  3)  and  (B.  5)  imply  the  familiar  result 

,2  2 

v(q)Il *  (?)-♦-»  ~* - 5 —  as  |  §  |  -»  00  . 

q  q  1 


In  a  similar  manner,  one  verifies  the  following  sum  rules  for 


A(E)  , 
rr 


CO 

’  dEA(E)  ,  =  6  ,  , 

J  rr  rr 


(B.  6) 


J  dEE  A(E)rr/=  Tr{e"P(H“^N"n)[ar  [H,a*,]]  } 

_  co  )  4- 


(B.  7 ) 


If  the  basis  jr)  is  that  which  diagonalizes  this  last  relation  becomes 

CO 

f  dEEA(E)  ,  =  e  6  , 

J  rr  r  rr 

-00 

+  )  [(rsjvlr's')  -  (rs  |  v|  s  'r  '  )~\  Tr{e  ^  ^a^  a  ,} 

t-j  /  s  s 

ss 
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e6  ,  +  Y  [  (r  s  |  v  |  r 's  ' )  -  (r  s  |  v  |  s  '  r  ' )]  •  —  V  S '  ,  (z  / )  e 

r  rr  L>  ,  11  11  p  L>,  s  s  L 

C  C  •  f 


v°+ 


-  «  .,-(1) 

e  6  /  +  G  ,  , 

r  rr  rr 


(B.  8) 


by  virtue  of  Eq.  (16).  Combining  (B.  6)  and  (B.  8),  we  have  for  |  z  |  — •  <=> 


S '  /(z)=“  -  6  .  +  ~  [e  6  ,  +G(11] 

rr  z  rr  Z  r  rr  rr 

2 


(B.9) 


Then  because 


G  ,(z)  =  (z-e  )6  [S '  (z)~  1  ]  ,  , 

rr  r  rr  rr 


it  is  easy  to  see  that 


(B.  10) 


(i) 


G  /( z)  -♦  Gv  ,  as  z  -♦  00 
rr  rr  1  1 
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